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Abstract 



<D . 
Ft . We study the evolution of both scalar and tensor cosmological perturbations in a 

Randall-Sundrum braneworld having an arbitrary expansion history. We adopt a four 

dimensional point of view where the degrees of freedom on the brane constitute an 

open quantum system coupled to an environment composed of the bulk gravitons. Due 

to the expansion of the universe, the brane degrees of freedom and the bulk degrees 

of freedom interact as they propagate forward in time. Brane excitations may decay 

through the emission of bulk gravitons which may escape to future infinity, leading to a 

r~| I sort of dissipation from the four dimensional point of view of an observer on the brane. 

Bulk gravitons may also be reflected off of the curved bulk and reabsorbed by the 

^SJ ■ brane, thereby transformed into quanta on the brane, leading to a sort of nonlocality 

^ i from the four dimensional point of view. The dissipation and the nonlocality are 

-^ I encoded into the retarded bulk propagator. We estimate the dissipation rates of the 

CN ■ bound state as well as of the matter degrees of freedom at different cosmological 

epochs and for different sources of matter on the brane. We use a near-brane limit of 



r^ I the bulk geometry for the study when purely nonlocal bulk effects are encountered. 

00 

1 Introduction 

rS \ Randall and Sundrum introduced a braneworld model with an infinite extra dimension [T] 

pi \ where our observable universe is a four dimensional hypersurface, the brane, embedded into 

a five dimensional Anti-de Sitter (AdS) bulk spacetime. The curved geometry oi AdS space 
induces an effective compactification in the extra dimension through its curvature radius 
i and thus prevents four dimensional (4d) gravity from leaking into the extra dimension. 
Consequently, standard 4d Einstein gravity is recovered at large distances on the brane. 
Corrections to the Newton square inverse force law arise at scales smaller than i [2] (in 
current experimental tests of Newton's law, £ < 0.1 mm [3]). Spatially homogeneous and 
isotropic cosmological background solutions in the Randall-Sundrum (RS) model have been 
found |4j. These solutions include an additional term quadratic in the energy-momentum 
tensor in the Friedmann equation, causing significant modifications of the standard cos- 
mology at high density. In the language of string theory, all the Standard Model fields 
are confined to the brane. Only the gravitons propagate in the entire bulk spacetime. It 
is hoped that the presence of an infinite extra dimension may be tested by cosmological 
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observations, such as the measurements of temperature and polarization anisotropies of the 
Cosmic Microwave Background (CMB). 

Cosmological perturbation theory about the cosmological background is an essential 
tool to probe the presence of extra dimensions. In braneworlds, cosmological perturbation 
theory is more complicated than in standard 4d cosmology because the brane breaks the 5d 
spatial homogeneity and isotropy of the bulk. Consequently, the linearized equations are no 
longer ordinary differential equations (ODE) in time but become partial differential equa- 
tions (PDE) in time and the extra dimension. Analytic solutions to cosmological perturba- 
tions have however been found for the maximally symmetric cases of a Minkowski brane or 
a pure de Sitter brane [I], [5]. In these special cases of higher symmetry the spectrum of 5d 
gravitons includes one discrete mode bound to the brane (if the bulk is Z2-symmetric) and 
a continuum of free modes. From the point of view of an observer on the brane, the bound 
state is interpreted as the 4d massless graviton and reproduces standard 4d Einstein gravity 
on the brane. The continuum is interpreted as 4d massive Kaluza-Klein (KK) gravitons 
whose interaction with the brane degrees of freedom is suppressed at low energy. In order 
to explore observational signatures for the presence of extra-dimensions, it is necessary to 
compute the evolution of cosmological perturbations in the more realistic case of a brane 
having a time-dependent stress-energy, i.e. a Friedmann-Robertson- Walker (FRW) brane, 
for instance with an inflation era followed by a radiation/matter dominated era. Several 
formalisms exist for evolving linearized perturbations in a braneworld scenario [7]. Here we 
use the formalism of Mukohyama [8]. The main difficulty is that, except for the case of an 
uniformly accelerated brane (flat brane or de Sitter brane) where analytical solutions have 
been obtained, the equations are generally not separable because of the complicated motion 
of the expanding brane in the bulk (a review on cosmological perturbations in expanding 
braneworlds is given in reference |12)). This problem has been investigated numerically 
and a number of interesting results have been obtained [T4l [T5l [T6l [TTl [T8l [T9l[20l[2n 122] . 
The braneworld perturbations problem is in principle entirely solvable numerically once 
the initial vacuum in AdS is known. However the specification of initial conditions in the 
bulk poses a problem of more fundamental nature. There is no unique physically motivated 
choice of initial conditions for an AdS bulk background [24l [25] (quite unlike the situation 
for dS space). Consequently the results of these numerical studies are subject to certain 
assumptions concerning the initial conditions. 

Five dimensional transverse and traceless metric fiuctuations describe the true bulk 
gravitational waves (5d gravitons) and contain five degrees of freedom (d.o.f). By "project- 
ing" the five d.o.f onto a fiat static brane, they split into tensor perturbations (the two d.o.f 
of the 4d graviton), vector perturbations (the two d.o.f called "graviphotons") and scalar 
perturbations (one d.o.f called "graviscalar"), in a S0{3) representation. This "projection" 
is somewhat ambiguous in the case of an arbitrary motion of the cosmological brane. In 
this case it is more convenient to use directly the 5*0(3) symmetry of the unperturbed 
spacetime to decompose the linear perturbations into tensor, vector and scalar perturba- 
tions. The advantage of this decomposition is that scalar, vector and tensor perturbations 
evolve independently at linear order. Tensor perturbations are transverse and traceless 
and are the easiest to deal with. They correspond to 4d gravitational waves and propagate 
freely into the bulk, independently of the presence of matter on the brane. By definition 
tensor perturbations are automatically gauge-invariant at linear order. Their wave function 
obeys an equation similar to a Klein-Gordon equation for a massless scalar field minimally 
coupled to gravity. The Israel junction conditions reduce merely to Neumann boundary 



conditions on the brane. Vector and scalar metric perturbations are more complicated to 
deal with because of their coupling to matter on the brane. Mukohyama has shown that, 
in the absence of matter in the bulk, scalar (as well as vector) perturbations in the "5d 
longitudinal" gauge can be derived from a single scalar master field, which obeys a five 
dimensional wave equation [8] (see also [9]). For vector perturbations, the boundary condi- 
tions for the master field are Dirichlet boundary conditions. In the more interesting case of 
scalar perturbations, the master field satisfies more complicated "nonlocal" boundary con- 
ditions on the branqll. Some authors have shown the connection between the master field 
formalism and the formalisms in different gauges ^0\. Although there has been progress in 
braneworld cosmological perturbation theory, there are a few quantitative predictions for 
the evolution of cosmological perturbations in expanding braneworlds. 

In this paper we study the evolution of cosmological perturbations of a FRW brane 
having an arbitrary motion in the AdS bulk. The goal of this work is to estimate the order 
of magnitude of the extra dimensional effects. Instead of looking for exact solutions to 
the perturbation equations, we are primarily interested in the role of the bulk gravitons in 
the brane-bulk interaction, following the ideas developed in the article of Binetruy, Bucher 
and Carvalho [26]. When the Hubble parameter on the brane changes with time, the ac- 
celeration of the brane in the bulk also changes and gravitons are emitted into the bulk. 
These may escape to future infinity or they may be reabsorbed by the brane because of 
refiections of the emitted gravitons off of the curved AdS bulk. From the point of view of 
an observer on the brane, these processes appear to generate dissipation and nonlocality 
[26] . Here these bulk induced effects are encoded into the retarded bulk propagator, which 
is systematically inserted into the effective brane propagator by doing a resummation of the 
bulk backreaction effects at all order in the brane-bulk coupling. The approach used in this 
work is a four dimensional perspective that regards the degrees of freedom localized on the 
brane universe as an open quantum system coupled to a large environment composed of the 
bulk gravitons. In the language of non-equilibrium quantum field theory, the bulk propa- 
gator plays the same role as a self-energy, dressing the bare fields composed of the discrete 
degrees of freedom localized on the brane. In this work we estimate the magnitude of the 
decay rate of the graviton bound state by studying the evolution of tensor perturbations 
in the frame of a cosmological brane with an arbitrary expansion history embedded in a 
Z2-symmetric AdS space. Some authors obtained analytical results for the order of magni- 
tude of the modified power spectrum for tensor perturbations, in the simplified braneworld 
model where the Hubble factor on the brane changes instantaneously [Tl]. Here we consider 
a more realistic model where the Hubble factor changes continuously and adiabatically in 
the sense that H <^ H^ . We also work out the dissipation of purely localized brane degrees 
of freedom, such as an adiabatic perfect fiuid or a slowly rolling scalar field on the brane, 
by studying the evolution of scalar cosmological perturbations and the coupling between 
metric and matter perturbations. We use Gaussian normal (GN) coordinates, where the 
extra dimensional coordinate measures the proper distance from the brane and the position 
of the brane in the bulk is fixed. GN coordinates have also the advantage of simplifying 
the form of the perturbed boundary conditions on the brane. The main drawback of GN 
coordinates is the presence of coordinate singularities in the bulk at a finite distance from 



^The term "nonlocal" for the boundary conditions has nothing to do with the nonlocality studied in this 
paper which is due to the information coming from the bulk. As it was stressed in [6], the terminology 
"nonlocal" for the boundary conditions is not appropriate, although it is frequently used in the literature, 
because they contain only a finite number of derivatives. 



the brane, despite the bulk space being regular and extendable beyond the singularity by 
choosing another set of coordinates. Because of the arbitrary motion of a FRW brane 
in AdS space, the metric components in GN coordinates have a complicated form which 
renders the equations of motion not separable. That is why we use a near-brane limit to 
perform the separation, actually only when purely nonlocal interactions arise between the 
bulk degrees of freedom and the brane degrees of freedom that are considered. This limit 
can be legitimized in the Randall-Sundrum model because the support of the wave function 
of the bound state is localized near the brane. We think that physics near the brane is 
suitable to describe dissipative effects in braneworlds. In order to focus our attention on 
the dissipative effects, we may approximate at high energy {H£ ^ 1) the inhomogeneity 
of the bulk responsible for backscattering of gravitons in the bulk and subsequent nonlocal 
effects on the brane. This approximation is discussed in section 16.11 

The paper is organized as follows. In section [2] we introduce the dissipative and nonlocal 
effects which appear in braneworld cosmology and we present the self-energy approach 
used in the paper to study these effects. In section [3] we briefly recall the cosmological 
perturbation theory in braneworlds. In section[4]we explore the scalar perturbations and the 
brane-bulk interaction for a slowly rolling scalar field on the brane. We apply our method 
to compute the local dissipation rate of the inflaton perturbation due to the interaction 
with the bulk. In section [5] we compute the Gaussian normal metric of AdS space with 
general Friedmann-Robertson- Walker (FRW) slices and present the near-brane limit of 
this geometry, which will be used as an approximation to evaluate the nonlocal corrections 
arising in the brane-bulk systems considered in the next sections. In section [6] we study 
the tensor perturbations and compute the dissipation rate of the graviton bound state. In 
section [7] we study the scalar perturbations for an adiabatic perfect fluid on the brane. We 
present our conclusions in section [8l 

2 Dissipation and nonlocality in the expanding braneworld 

Cosmological perturbations theory for the Randall-Sundrum model has been solved ana- 
lytically for special brane universes where the Hubble factor on the brane H is constant. In 
these cases the four dimensional universe is static {H = 0) or follows a uniform expansion 
{H > 0). In the flrst case the Minkowski brane is at rest with respect to the five dimen- 
sional Anti-de Sitter (AdS) bulk, whereas in the second case the pure de Sitter (dS) brane 
follows a uniformly accelerating trajectory in the AdS bulk. The equations for the pertur- 
bations may be reduced to a Schrodinger-like equation. The presence of the brane within 
the bulk creates a (5-function potential well, while the curvature of the AdS bulk creates 
a decreasing barrier potential (Fig. [T]). The resulting volcano-like potential gives rise to a 
single brane bound mode, the zero mass mode m = 0, in the spectrum of the linearized 
Einstein equations, surrounded by a continuum of massive Kaluza-Klein free modes, with 
m > 3H/2 for a dS brane [5]. The zero mode represents the standard massless graviton 
and reproduces the four dimensional gravity on the brane. The amplitude of the massive 
gravitons is suppressed near the brane at low energy due to the barrier potential. 

Because of the presence of matter in the universe, the cosmological expansion history 
is no longer uniform. The Hubble parameter changes with time, inducing a change in the 
shape of the potential. Consequently, the brane degrees of freedom and the bulk degrees 
of freedom interact thus generating transitions between the modes. 
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Figure 1: Left: Volcano potential for a Minkowski braneworld. The height of the potential 
is 0{i~'^) and the potential decreases as 1/z^ if z is the extra dimension. The zero mode graviton 
bound to the brane is surrounded by a continuum of massive gravitons m > 0. Right: Volcano 
potential for a de Sitter braneworld. The height of the potential is 0{H^ + £~^) and the 
potential decreases as 1/sinh (z).The zero mode graviton bound to the brane is surrounded by a 
continuum of massive gravitons with m > 'iH/2 characterized by a gap. 



The brane-bulk system is a Hamiltonian quantum system, which is necessarily conser- 
vative because the phase space density must be preserved under time evolution. In this 
sense, the brane-bulk system is intrinsically non-dissipative, and the appearance of dissipa- 
tion can arise only as the result of "coarse graining". The interaction between the quanta 
due to the expansion of the brane may be characterized by a Bogoliubov transformation 
[S matrix for a linear system) relating the modes of positive and negative frequency for 
the initial "in" vacuum and the resulting "out" vacuum. However, from the point of view 
of an observer confined to the brane, the four dimensional universe as a submanifold is 
an open quantum system. Consequently, the mixing between the modes localized on the 
brane and the delocalized modes in the bulk creates the appearance of dissipation to a 
four dimensional observer unable to access the bulk modes. The observer on the brane is 
sensitive only to a part of the complete AdS vacuum, composed of only the discrete modes 
localized on the brane. We call this part "brane vacuum". The Hilbert space is truncated 
because of dimensional reduction, and the "bulk vacuum", composed of the continuous bulk 
modes, holds the missing information. When we observe the cosmological perturbations 
today, we measure expectation values of observables quadratic in the creation and annihi- 
lation operators for the modes localized on the brane today, namely ahrane,out and a 
The S matrix expresses the "out" operators as linear combinations of abran 
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on one hand, and of abuik,ini^) ^^^ CLluikini^) ^^ ^^^ other. A useful parameterization 
of this transformation has been proposed in |26]. The authors require that Abrane,in and 
^buik,in be entirely on the brane and in the bulk, respectively, and be normalized such that 
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1. Then abrane,out may be expressed in terms of 



these according to one of the three following possibilities: either 
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where < ^ < 7r/2; or 



Q'brane,out — COSll UAbrane^in ~r Smil U-n.^ 



bulkAn 



(1) 

(2) 



where < n < +cxd; or 

abrane,out = sinh uAI^^^^^^^ + COsh uAbulk,in (3) 

where < n < +00. AbraneM may be constructed entirely as a linear combination oiabrane,in 
and flferanein' ^^'^ likcwisc Abuik,in may be constructed entirely as a linear combination of 
the abuik,in{k) and the aluikinik)^ ^ being the index for the continuum modes. We observe 
in ([31) that the bulk initial state may have a very important, or even dominant, role in 
determining what we see on the brane today. 

However it has been numerically observed in [H] that the evolution of tensor pertur- 
bations on the brane today does not show any significant dependence on the bulk initial 
state when a rfS'-invariant initial vacuum is specified in the AdS bulk, that is defined in 
a (iS'-slicing frame. However we might suspect that the initial amplitude of bulk modes 
that are initially located outside the bulk Cauchy horizon of the (iS'-slicing will not be 
suppressed when they causally affect the decelerating brane in the future. The relative 
contribution of the brane initial conditions and the bulk initial conditions in determinining 
what we see on the brane today should depend on the basis decomposition of the m-vacuum 
as described in ([3]). The main problem is that there are no natural initial conditions for 
AdS [in [25]: in standard inflation the backgound geometry is dS so the timelike geodesies 
diverge forward in time losing causallity such that initial irregularities become swept out, 
dS space thus exhibits natural initial conditions for generating homogeneity and isotropy. 
Whereas in Randall-Sundrum branewords the background geometry is AdS so that the 
timelike geodesies first diverge and then refocus forward in time, thus conserving the initial 
amplitude of metric perturbations. 

The brane-bulk interaction may be summarized by the following fundamental processes 
from the four dimensional point of view illustrated in Fig. [3 An initial vacuum state may 
be completely characterized by specifying the quantum state of the incoming gravitons on 
the past Cauchy horizon in the bulk if(_) and of the degrees of freedom on the brane at 
the intersection of the brane with -ff(-). Subsequently, the bulk and the brane degrees of 
freedom interact as they propagate forward in time. Bulk gravitons may be absorbed and 
transformed into quanta on the brane. Similarly, brane excitations may decay through the 
emission of bulk gravitons. These may either escape to future infinity, leading to dissipation 
from the four dimensional point of view, or may be reabsorbed by the brane due to the 
bulk curvature, leading to nonlocality from the four dimensional point of view. 

We may consider the brane-bulk interaction from the four dimensional perspective by 
regarding the brane degrees of freedom as an open quantum system coupled to an environ- 
ment with a large number of degrees of freedom. Schematically, we may block diagonalize 
by Fourier transforming in the three transverse spatial dimensions. For an expanding 
braneworld with a perfect fluid or scalar field on the brane, the interaction between the 
scalar brane degree of freedom q{t) and a scalar bulk degree of freedom u{t, x) may be rep- 
resented by a system of coupled equations. Under the approximations presented in detail 
later in the paper, this system of coupled equation may be written according to the general 
form 

[-dta(t,x)dt + d^h{t,x)d^ + c{t,x)]u{t,x) = 0, 

[d, + \{t)]uit,x)\,=o = Piidt)q{t), 

[d^+u^,{t)]q{t) = P2{dMt,x)U=o, (4) 



H(+) 



bulk 



H(-) 



timelike 
brane 



Figure 2: Fundamental processes for the brane-bulk interaction from a four dimen- 
sional perspective. Prom the left to the right: absorption, dissipation, nonlocality. The nonlo- 
cality arises from the diffraction of gravitons into the curved AdS bulk. 

where Pi and P2 are polynomials in dt of degree one and degree two respectively, with time 
dependent coefficients. The first equation is the bulk equation of motion, the second the 
boundary condition on the brane, and the third the equation of motion for the brane degree 
of freedom. The brane and the bulk degrees of freedom are related according to 
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dt'Gtrane{t,t')P2{dt')u{t',0), 
dt'Ghulk{t,t',X,0)Pi{dt,)q{t') 



(5) 



where Gbrane = [d'i ~^ ^o{^)] is the barcl retarded brane propagator and Gbuik is the 
Neumann form of the retarded bulk propagator. Due to the interaction, the effective 
retarded brane propagator may be obtained by resumming the infinite geometric series at 
all orders of the bulk backreaction 
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(6) 



where the bulk propagator links two points both lying on the brane. The brane degree of 
freedom is governed by the integro-diff'erential equation 






brane 



dlq + ul{t)q + / dt'Kbuikit - t')q{t') = 0. 



Here the bulk kernel is given by 

Kbulk(t,t ) 



-P2{dt)Gbuik{t,t')Pi{dt> 



(7) 



(8) 



and plays the role of a self-energy, dressing the bare field on the brane. Both the local 
dissipative effects and the nonlocal effects due to the interaction with the extra dimension 



^By "bare" we mean the propagator of the free brane fields excluding the interaction with the bulk. 



are contained in this self-energy kernel. The bare retarded brane propagator Gbrane = 
[d^ + u!Q{t)] is the propagator of an oscillator localized on the brane and accounts only 
for the transitions between the four dimensional modes due to the expansion of the universe. 
We consider as an illustration the following couplings 

Pi{dt) = -li{t)du P2{dt) = -i2{t)dt. (9) 

Since the Neumann retarded bulk propagator on the brane can be written in the form 

Gbuik{t. t', x = 0,x' = 0)= 9{t - t')G{t, t', 0, 0), (10) 

the bulk kernel may be decomposed into singular and regular parts 

K,Mt,t') = K-^{t,t') + KZUt,t') (11) 

where 

KZ7it,t') = d{t-t'h,{th2{t)G{t,tAm 

Kuikit.t') = e{t-t')^2{mG{t,t',o,o)j,{t')d,. (12) 

The singular part is responsible for the local dissipative processes and the regular part 
describes subsequent nonlocal processes, such as reflections from the inhomogeneous curved 
bulk. The infinite summation 

^^brane ^^brane ~r '^brane \ ^^bulk j ^-'brane ~T ^^brane \ ^^hulk j ^^brane \ -'■^ bulk J ^^brane ~r ..., 

is equivalent to adding a local dissipation term r(t) = ■yi(t)'j2(t)G(t, t, 0, 0) to the equation 
of motion for the brane degree of freedom. The infinite summation 

'^brane ^brane ~r ^brane \ ^bulk) ^brane 

+Gbrane {~-i^bulk) Gbrane \~i'^bulk) Gbrane + ••• (14) 

adds the nonlocality to the equation of motion such that the brane degree of freedom 
propagates as 

G^\ne°'i = m+mm+^iiMt)+^2{t) f rft'G(t, t',o, 0)71 (t')g(o = 0.(15) 

Since a local singular part does not depend on the curvature of the background, the term 
G*(t,t, 0,0) which appears in the local dissipation rate may be replaced by the Neumann 
form of the Minkowski propagator. We recall that the Minkowski propagator is 

GM^nk{t,t',X,x') = 6 ((t - tf - {x - x')^) Jq (pV{t " ^0' ' {x - X')A , (16) 



where p is the momentum in the three transverse spatial dimensions. Then G(t,t, 0,0) = 
GMink{t,t, 0,0) = 1 and the local dissipation rate depends only on the couplings 

m = 7i(t)72(t). (17) 



For certain derivative couplings, additional local terms may appear leading for instance 
to a phase shift. There might also be a nonlocal dissipation term when there is no derivative 
in the couplings. This formalism with Green functions is useful to discriminate transitions 
among four dimensional brane modes due to the expansion of the universe from transitions 
between the brane and the bulk modes due to the presence of an extra dimension. Moreover 
this formalism allows us to distinguish the local dissipation processes from the nonlocal 
processes. 

The purpose of this work is to estimate the dissipation rates of certain degrees of freedom 
confined to the brane or localized near the brane from an arbitrary expansion history of a 
brane having no special symmetries. 

3 Braneworld cosmological perturbations 

In this section we recall briefly the cosmological perturbations equations in braneworlds and 
the formalism of Mukohyama frequently used in the literature [6l [U [9l [IHl [131 EH [IE! [22]. 
The metric perturbations hf^u about the general bulk background metric g^u, in Gaussian 
normal coordinates, given by 

ds^ = g^.^dx^dx'^ = d^^-N\^,T)dT^ + A\^,T)dxl (18) 

describe the five dimensional bulk gravitons. The evolution of these bulk perturbations 
(where we assume /i^,y <C g^u) will depend on the coupling to the matter perturbations 
ST^i, localized on the brane at ^ = through the linearized Einstein equations 

SG^, (V) = ^^^ST^u (19) 

from which linearized Israel junction conditions on the brane can be derived. The coupling 
is related to the five dimensional Planck mass k^ = Stt/M^; = SttGat^, where Gn is the four 
dimensional gravitational coupling constant and £ the curvature radius of the bulk space 
in the Randall-Sundrum model. The bulk metric perturbations can be decomposed into 
scalar, vector, and tensor perturbations. 

3.1 Tensor perturbations 

The gravitational waves away from the brane propagate freely in the bulk and independently 
of the presence of matter on the brane. They are described by transverse and traceless 
tensor perturbations Eij of the metric 

rfg2 = (^^^ + /^^^) dx^dx^ = di'^ - A^2(e, T)dT'^ + A^{i, t) {5ij + Eij) dx'dx^ . (20) 

Tensor perturbations are gauge invariant at linear order and they do not couple to matter on 
the brane. The tensor perturbations are exactly described by a minimally coupled massless 
scalar field \l/(t, ^;p), Fourier transformed in the three transverse spatial dimensions, as 

E,,(t,e,x) = Jdp^it,^;p)e'P-^ei0e„ (21) 

which propagates in the bulk background according to the Klein-Gordon equation 
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^^ = 0. (22) 



The Israel junction conditions reduce to a Neumann boundary condition on the brane 



a^ 



?=o 



0. 



(23) 



This boundary condition is homogeneous under the assumption that there are no anisotropic 
stresses on the brane. 

3.2 Scalar perturbations 

The scalar perturbations of the bulk background metric may be simplified by choosing a 
particular gauge knowned as "5d longitudinal gauge". In this gauge there are four gauge- 
invariant variables and the line element is 



ds 



■^2 



(1 + 2A^i:) d^^ - N\^, t) (1 + 2 A) dr^ + N{^, r)A^d^dt 
+A\^, r) (1 + 27^) 6ijdx'dx\ 



(24) 



These five dimensional scalar metric perturbations couple to the four dimensional scalar 
matter perturbations on the brane, which are defined in the same gauge on the brane b}0 
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/lU 
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flU 



-(P + Sp) 6qj 

-A~Hq^, (P + 5P) 6^j 



Here the greek indices label the four dimensions on the brane and the latin indices label the 
three transverse spatial dimensions on the brane. Throughout the paper, p and P denote 
respectively the energy density and the pressure of the effective matter on the brane as 
opposed to the real matter on the brane labeled by a subscript M and related to the 
previous one by p = pm + cr, P = Pm — o", where a is the brane tension, fine-tuned in 
the Randall-Sundrum model so that a = Q/{n^tj [I]. Mukohyama [8] (see also [9]) has 
shown that in the absence of bulk matter perturbations, the scalar metric perturbations 
in 5d longitudinal gauge can all be generated from a single scalar "master" fielcu f2(t,^) 
according to 
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(25) 



The master variable fl satisfies the five dimensional wave equation 
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^We have neglected the scalar anisotropic stress which do not arise when considering the perturbations 
of a perfect fluid or a scalar field. 



■^It should be noted that [n] = {length^. 
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The Israel junction conditions reduce to the following "nonlocal" boundary conditions on 
the brane [T3l 
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(27) 



4 Scalar perturbations: slow-roll infiaton on the brane 

In this section we work out scalar metric and matter perturbations for a scalar field on 
the brane with a slow-roll potential, V{(f>), such that the induced geometry on the brane 
is quasi-de Sitter. During slow-roll inflation the expansion of the universe is adiabatic in 
the sense that H <^ H'^. The infiaton scalar field, 0(r), may be characterized by an energy 
density Pm and a pressure Pm 



1 > 



Pm 
Pm 



Vi<l>). 



(28) 



We may combine the equations of Mukohyama presented in section 13.21 for scalar per- 
turbations to obtain simple coupled equations of the brane-bulk system. Following the 
calculations done in [lH [22], we introduce the Mukhanov-Sasaki variable as the brane 
degree of freedom: 



Q 



^7^„ 



(29) 



where S(f> is the perturbation of the infiaton and TZb is one of the scalar metric perturbations, 
namely the curvature perturbation (1251 ). projected on the brane. The equations of the 
brane-bulk system in the general Gaussian normal coordinate system ( ITSl ) are then given 
by dUEg 
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n- 3 



A N 

A^N 



Q + 3HQ 



p^ 
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Q" 



'^A iv'^^ 
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Hin'-^n)^^-^ 



V"{<P) Q 




H Hl^ (H 
H H A, \H 



where p is the transverse momentum and the source J(f2) is given by 



2 i2 

K ad) 



H 



-Q 



J(fi)l€=o, (30) 
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f]' 



ij \ p-^n 



2n'^ / 9a5 
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n^6a^ 



Here the notations are a^r) 
equation l \102f \ we also have 



nQ'-HQ--Q + j-^Q 



(31) 



A{t,^ = 0) and 7i = {A'/A)\^=o. From the Friedmann 



K^072 



-iH/Vi+¥iP. 



(32) 



The brane-bulk equations (l30l ). (l3T1 ) may be simpUfied for slow-roll inflation as follows: 
we neglect all the adiabatic corrections to the expansion, like H/H"^ terms, except for the 
terms involved in the coupling between Q and Q. According to these approximations and 
after rescaling the master field according to AQ ^ Q, we argue that the system of coupled 
equations reduces to 



Qir') 



^braneVi '^ 



0, 

Pi{dr)Q{r), 

P2(5.)fi(r,0ls=o, 



(33) 



where Gbuik is the retarded Green function of exact bulk equation of motion, 
Gbrane = [9'^ + 3H{T)dr + p^ / tt^] is the bare Green function of the inflaton perturbation 
onthebrane, Pi(a,) = -{K^a^^/p^)dr d.TvdP2{dr) = -{p^/{6a^H)) [d^ + H{T)dr + p^a^]. 
The first equation in ([33l) is the bulk wave equation. The presence of time derivatives of the 
fields in the brane-bulk couplings induces local dissipation of the scalar field on the brane. 
The effective brane propagator with the interaction with the bulk taken into account is 
obtained as explained in section [2] by resumming the infinite geometric series 



'~^ brane 
'^ brane 



G. 



brane 



p^(j) 



^bulk 



^^bulk "I o^bulk 

a'' 



^-1 
brane 



Ga'^H ['^bulk 



^^bulk 



"I" a'^^bulk 



-Dr 



K^a (f) 

p2 



Dr 



G. 



brane 



(34) 



N 
bulk 



where G^^^^^ = [D^ + SHD^. + p^/a^] describes the bare propagation on the brane. G 
is the Neumann form of the bulk retarded propagator projected on the brane and has the 
form 



GZik{r,r', ^ = 0,e = 0) 



r-r')G^ir,T', 0,0). 



(35) 



The time derivatives of the bulk propagator ([35l) appearing in the effective brane propagator 
( l34l) lead to singular and regular terms. This suggests that the dressed propagation of the 
infiaton contains local and nonlocal adiabatic corrections 
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G^^(r,r, 0,0) 
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3^ '^ / zjr^Nt^ ^ n n\ , r^N , 
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HG^{t,t,0,0) + G'^t,t,0,0) 
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Q + (nonlocal term) = 0. 



(36) 



Here the nonlocal term depends on the AdS curvature i and the brane intrinsic curvature 
H, and is given by the time derivatives of the regular part G^ of the AdS bulk Green 
function G^^^. of the bulk wave equation: 



2a2 



(nonlocal term) 






ds 



G''{s) + HG''{s) + ^G''{s) 



Q{r-s). (37) 



The local terms do not depend on the curvature, which means that the regular part 
G^ [t, t, 0, 0) of the bulk Green function is equal to the Minkowski propagator at the origin 
in ( l36l ). Since the Neumann form of the Minkowski propagator Gj^j„^(r, r',^,^') is given 
by 






Mink 



(r, r', e, O = GM,^^k{T, r', ^, O + G,unk{r, r', ^, -^ 



(38) 



where 



GM.n.(r,r',e,n = {{r - r'Y - (^ - e?) -Jo (p^/F^W^If^W) , (39) 



one has G^(r,r, 0,0) = Gj^/i„fc(r, r, 0, 0) = 1 and G^(r,r, 0,0) = Gj^/,„fc(r, r, 0, 0) = so 
that the effective equation of motion for the inflaton reduces to 
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Q 



6H 



Q + '^Q + (nonlocal term) = 0. (40) 



We may renormalize the kinetic term to one by dividing this equation by the coefficient in 
front of the kinetic term to obtain 



Q+[3H + r(r)] Q 



p" 



A(r) Q + (nonlocal term) + 0(kV ) = 0. (41) 



From (14T|) we observe that the interaction of the inflaton with the bulk gravitons lead to 
local dissipation into the extra dimension through the local friction term appearing in the 
effective equation as the first time derivative V{t)Q. There is in addition a phase shift 
given by K{t)Q. We find that the local dissipation rate of the infiaton perturbation due to 
the extra dimension is: 



2 J,2 



Tir) 



K^(j) 



0{l) 
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He 



H^VTTlPP 
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and the local phase shift is 



A{i 



(9(1) 



H 



Hi 



p^ 



H'^^/l + HH^a?' 
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(42) 



(43) 



We observe that the dissipation term TQ dominates the phase correction AQ at superhori- 
zon scales and may be comparable at high-energy {Hi ^ 1) to the slow-roll corrections to 
standard inflation. If inflation takes place in a regime where the bulk curvature radius is 
much larger than the Hubble radius {Hi ^ 1), then the local dissipation rate behaves as 

r(r) ^'^>^ 0{1)^^H (44) 

In a quasi- four dimensional regime {Hi <^ 1) the local dissipation rate behaves as 

r(r) ^^<^ 0{l)^^{Hi)H (45) 

The dissipation rate is suppressed in a quasi- four dimensional regime by the factor {Hi). 
The scalar field on the brane dissipates linearly with the slow-roll parameter at any scale. 
As we will see for tensor perturbations, the brane-bulk interaction between the graviton 
bound state (zero mode) and the continuum of bulk gravitons will be purely nonlocal since 
the brane degree of freedom here is not confined but localized near the brane. In case of 
purely nonlocal interaction we need to consider the bulk background geometry, which is 
the object of the next section. 

5 Background metric 

A flat (3+l)-dimensional homogeneous and isotropic Friedman-Robertson- Walker universe 
with an arbitrary expansion history, characterized by the scale factor a{T) as a function of 
the proper time r, can be embedded in a wedge of AdS^ of curvature radius £, described 
in terms of the static bulk coordinates (Poincare coordinates) 

ds^ = f -dt^ + dfl + d^\ (46) 

by means of the following explicit embedding 



dr 



z,{t) = e-J -^-^'^^ U{t) = J -^Ji'zAr) + z!{r), (47) 

where H{t) is the Hubble factor on the brane. We may explicitly construct Gaussian normal 
coordinates by computing the geodesic curve in the zt-pla.ne normal to the braneworld 
trajectory at proper time r, traveling at a proper distance ^ away from the brane. We 
thus obtain the following mapping from the Gaussian normal coordinates to the Poincare 
coordinates: 

p-JdT H{r) 
z{i,T) - 



cosh(e/£) - Vl + ^2^2sinh(e/£)' 

iHe-i'^^ ^^^Uin\i{i/i) 



t{^,T) = /'^[e-/'^-^MVl + £2^2 



from which follows the line element 



COsh{^/i) - y/l+im^{T)smh{^/i) 



,(48) 
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ds' 



VlT¥IPsmh{^/i) - cosh{^/i) + 



'H 



VT+WlP 



smh{^/i) dr^ 



+de. 



(49) 



In these coordinates the brane is stationary with respect to the bulk, and C, measures the 
proper distance from the brane. 

Despite being convenient for describing a brane with an arbitrary expansion history, 
Gaussian normal coordinates suffer from a number of drawbacks. In particular this coor- 
dinate description can break down whenever the spatial geodesies normal to the brane de- 
velop caustics, either by focusing in time at ^h = ^coth~^ f \/l + PH"^ + i'^H/y/l + PH^) , 

or in the transverse spatial dimensions at ^s = ^coth^^ (Vl + (?H'^)- Even though these 
singularities may give the illusion of a horizon, the bulk space is regular there and can be 
extended beyond by using another set of coordinates. Let us consider the initial value prob- 
lem in these coordinates by looking at the Carter-Penrose diagram for Randall-Sundrum 
universes (Fig. [3]). We specify initial data on a surface of constant Gaussian normal time, 
in the bulk as well as on the brane. Suppose that we limit our ambition to predict what 
happens on the brane in the future. The evolution of a dust-dominated braneworld uni- 
verse, which was initially inflating in the past, may be causally affected in the future by 
unknown information coming from outside the initial past Cauchy horizon, because the size 
of the bulk horizon has increased during the deceleration of the expansion. 



z = 



z = infinite z = 



z = infinite z = 



z = infinite 






Figure 3: Carter-Penrose diagrams for various braneworld cosmologies embedded into 

AdS. The heavy line indicates the trajectory of the timelike brane. The horizontal direction 
represents the extra dimension. The past Cauchy horizon is noted H{—). From the left to the 
right: Minkowski braneworld, de Sitter braneworld, Friedmann-Robertson- Walker braneworld. 
Panel on the right describes an universe whose expansion initially inflates and then slows down. 

The complicated arbitrary motion of the brane into the bulk breaks the time-translation 
symmetry of the bulk, such that the metric components in eqn. (|49|) exhibit a non- 
separable form. The evolution of the bulk perturbations in the background metric (j49l) 
is not amenable to analytic methods without some form of approximation. Consequently, 
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we will modify the metric retaining only its most important features and assume that the 
bulk solution is quasi-separable. From the point of view of an observer on the brane, most 
of the action takes place within a small number of AdS decay length (or apparent decay 
length) of the brane, where virtually all the bulk four-volume is concentrated. The bulk 
modes bound to the brane have most of their support localized there, therefore a mediocre 
approximation of the metric far from the brane is only likely to provide a poor approx- 
imation of the tail of the wave function, where there is almost no probability. We may 
also hope that any quanta that escapes from the brane, due to non-adiabatic effects in the 
expansion history, do not come back through reflection or diffraction. Consequently, in 
this case an observer on the brane will care little about how exactly such quanta escape, 
which will depend on the metric away from the brane. In the WKB approximation, most 
of the quanta become classical within a short distance from the brane. In that sense we 
may approximate the background geometry ([49l) by the geometry near the brane {^ <^ i), 
by means of the following line element 

ds"" ^ de - e-2"^M^rfr2 + e-2"^«^a2(r)dx2, (50) 

where we have Z2-symmetrized about ^ and where the time dependent warp factors and 
the scale factor are respectively given by 



ai(r) = Vl + PH^ + 



Vr+72^ 



a2(r) = Vl + PH\ 
a{T) = e^^(")''". (51) 

In this approximation (l50l) the coordinate singularities are also removed to infinity. 

6 Tensor perturbations in the approximate background 

We now study the simplest case of tensor perturbations, which evolve independently of the 
matter content on the brane. We use the approximate background geometry (1501 ). which 
is accurate near the brane. Each polarization of the tensor perturbations is described by a 
minimally coupled massless scalar field \E', as discussed in section [3l We take the Fourier 
transform in the three transverse spatial directions and separately evolve each Fourier mode. 

6.1 The plateau potential 

The equation of motion of the massless scalar field \E' in the background metric (1501) is 

^)kl)9. + e--'¥(a|-»9nK)(^ + ^)a5) 

(52) 

where p is the momentum in the three transverse dimensions. From this equation in the 
interval —00 < ^ < +00, we may extract the homogeneous Neumann boundary condition 
( l23l) at .^ = for the even modes. After rescaling the scalar field as 

$(r,0 = e-("^(")+=^"2("))ia(r)i^(r,0, (53) 
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the equation of motion becomes 
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-i^i(^-^j + o(iei,ier) 



$ = 0. 




(54) 



The equation of motion is not separable yet. However, because the bound state is localized 
near the brane, an adequate approximation to the evolution of the bound state may be 
obtained by retaining only the leading behavior of the coefficients in ^ about ^ = 0, thus 
simplifying the equation. Doubtless this is a poor approximation for the tails of the bound 
state but we argue that the tails contribute negligibly. Consequently we set 
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:H' 



reducing eqn. (J54l) to 

with the effective plateau potential 

Vi^,T) = -VoiT)5iO 
where 
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p^ 



2 a\T) 
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(55) 



(56) 
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ViTPiP 



(57) 



(58) 



The corresponding boundary condition is obtained by integrating the Z2-symmetric equa- 
tion (l56l) . Equivalently, we may impose the boundary condition separately 



d,+ 



Voir) 



$ 



?=o+ 







(59) 



at .^ = and restrict the domain of (1561 ) to ,^ > 0, thus removing the (^-function. Under this 
approximation, the crater of the usual volcano potential [T],[5] is faithfully rendered but the 
landscape around the summit remains at constant elevation, as indicated in Fig. [4] right. 

When the brane has a pure de Sitter geometry, the linearized Einstein equations in 
AdS for the tensor perturbations may be reduced to a "Schrodinger-like" equation for a 
scalar field in the conformal bulk coordinates once the field has been rescaled. The exact 
expression of the effective potential has been calculated in [5]. The potential has the 
usual volcano shape with a decreasing barrier potential resulting from the curvature of the 
AdS bulk (Fig. [4] left). The height of the decreasing barrier potential is 0(£~^) and the 
energy gap between the zero mode and the continuous modes is 0(if^).Thus the plateau 
potential approximation we presented previously consists in neglecting the inhomogeneity 
of the bulk and consequently the diffraction of the bulk gravitons by the curved bulk. 
This approximation is legitimized in the cosmological regime where the curvature radius 
of the bulk is much larger than the Hubble horizon, H£ ^ 1. This is because in this 
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Figure 4: Approximation of the volcano potential by the plateau potential. 

regime the energy range of the barrier potential is of order £~^ and is much smaller than 
the Hubble energy scale H"^ of the gap. In this regime the volcano potential looks like 
the plateau potential (Fig. [4]), in which case the bulk inhomogeneity may effectively be 
neglected. In the regime where the curvature radius of the extra dimension is much smaller 
than the Hubble horizon, HI < 1, our approximation fails because the energy range of the 
inhomogeneity of the bulk is now non negligible. Thus in that regime the plateau potential 
would only give some rough lower and upper bounds of the true volcano potential. 

6.2 Decay of the bound state for an inflating brane 

The potential in (l57l) has exactly one bound zero mode and a continuum of free modes 
whose energy eigenvalues begin above the plateau. We expand the wave function in normal 
modes 



'0 



+ 00 



(60) 



where b labels the bound mode and k labels the continuum modes. The mode functions 

(pn{^]T) {n = b,k) satisfy 



-dl - Vo{t)S{0 + 



V,(t? 



0n(^;r) = m^(r)0„(^;r). 



The Z2-symmetric spectrum of (IBT! ) consists of a single bound state with niKr) = 0, 



.Mil 



0,(e; r) = iV,(r)e-^l«l, N^ir) = ^/W^ 
and a continuum of free bulk states k with m\{T) = k"^ + Vq{t)/A^ 



0fc(e;r) = iVfc(r) 
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(61) 



(62) 



(63) 



A gap of Vo(r)^/4 separates the bound state from the continuum. 

We consider the universe during an inflationary epoch, when the geometry of the brane 
is quasi-de Sitter and the expansion is adiabatic in the sense that H -C H^. In this regime 
the brane-bulk interaction may be computed perturbatively. The only tensor degree of 
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freedom on the brane is the discrete bound mode. We may compute the bound state decay 
rate due to the inflating motion of the brane into the bulk. Under the above approximation, 
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(64) 



The equation for the time dependent mode expansion coefficients Cnir) {n = b,k) in the 
approximate background described by eqn. (i50l) is 



Cn + 



a'^[T) 4 
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d4>n 



+ (higher order terms) 
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(65) 



where the higher order terms include two time derivatives acting on the potential (e.g. 
H,H^). The terms (0n|c^vb0n) vanish since the eigenmodes are real. The terms {4>k\dvo'Pk') 
connecting two continuum states are also zero because of orthogonalitja. The only nonzero 
matrix elements connect the bound mode and a continuum mode, 
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(67) 
(68) 



coupling the bound state to the continuum states with k > 0. In the adiabatic approxima- 
tion the first order coupling factor 7fc(r) and the frequencies Qb, ^k are 
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(69) 



The time dependent coupling factor in (l67l) leads to transitions between the bound mode 
and the continuous modes, because the zeroth order bare modes act as a source for the 
modes at next order. As soon as the acceleration of the brane changes in the bulk, through 
the derivative of the Hubble factor, H{t) (or equivalently Vo(r)), the in mode acts as a 
source through the coupling factor and generates different modes at first order. Conse- 
quently transitions from the bound mode to the continuum occur at first order, and the 
first order continuous modes act then themselves as a source for the bound mode, as a back 



^The orthogonality is the result of the simplified shape of the plateau potential. 
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scattering effect at second order. Tlie brane bound mode and tfie bulk continuous modes 
are related according to 



Ck[T) = I rfr'GL,(r,r')(-7.(r'))c,(r'), 

Cb{T) = I dT'Gi^ane{r,r') f dk^k{T')Ck{T'), (70) 

J —OQ J 

where G'^^ij.{t, t') and Gbrane{'T,'T') are the bare retarded Green functions satisfying 

{d'^+nl{T))Gbrane{r,T') = 5 {t - t') . (71) 

We may express the retarded bulk Green function of a mode k in the WKB approximation 

sin [J^, Qk{r)dT) 



GLkir^r') ^ 9{r-r') 



v/n^M^4py ■ 



(72) 



This WKB approximation is reasonable because that bulk eigenmodes oscillate at all times. 
There is no turning point because Vl\ ?» p^ /a? — 9H'^/4 + k"^ + Vq /A > always holds. The 
condition Clk <^ ^l is verified for both subhorizon modes, p ^ aH, and superhorizon 
modes, p <^ aH. The effective brane propagator with the interaction with the bulk modes 
of freedom taken into account is obtained by summing the infinite geometric series 
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(73) 



Gbrane + Jdk^kDrGl^.^kDr' 

The four dimensional interaction between the brane modes due to the expansion of the 
universe is contained in G^^^g, whereas the interaction between the bound mode and the 
bulk continuum is contained in the bulk kernel 

K{t,t') = [ dk^k{T)drGlikir,r')^k{r'), 

^ ^"f]fc(r)dr^7fc(r'), 



eir - r') / dk-fkir) 



fitfr) 



cos 



(74) 



such that the brane bound mode is governed by the integro-differential equation 






brane ^ 



d^Mr) + ^l{r)cb{T)+ f dT'K{T,T')dr'C,{T') = 



(75) 



where the bulk interaction kernel K{t, t') dresses the bare field of the bound state. The 
possible dissipation of the bound state appears to be nonlocal as viewed by an observer on 
the brane. This feature arises because the bound state is not a purely local but a localized 
brane degree of freedom with a typical extent in the bulk given by iatt = {Vq/2)~^. The 
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existence of the bound state is intrinsically dependent on the bulk curvature i~^. This 
condition is incompatible with a local interaction with the bulk. In the limit Hi ^ 1 one 
has Vq ~ AH, Vq = AH . The bulk interaction kernel (l74l ) reduces in the limit Hi ^ 1 to 
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at subhorizon scales, because the integral over k is suppressed ioi k ^ H in that case. 
We used that H, H are fairly constant in time in the adiabatic approximation. The 
superhorizon kernel ( f76l ) is plotted in Fig. El We observe that the kernel is nonlocal, 




Figure 5: Plot of the bulk interaction kernel K{s) with respect to time s, at super- 
horizon scales p « aH and in the limit Hi ;» 1. Here H = 1. 

typically varying with the Hubble time scale H~^. Thus the bulk interaction contains a 
memory. Consequently a local approximation would be possible only if the bound mode 
varied typically with a time scale larger than the Hubble time. 

During inflation, the Hubble radius of the universe remains fairly constant, whereas the 
physical wavelength of a mode increases exponentially due to the expansion. Early during 
inflation a mode of a given wave number is initially subhorizon so that p ^ a{T)H . As the 
universe expands, the mode crosses the horizon to become superhorizon. Whereas modes 
oscillate at subhorizon scales, they become frozen in after they exit the Hubble horizon. An 
observer is sensitive to the superhorizon modes on the brane which are the relevant modes 
to present cosmological observations and constitute the "initial conditions" of the universe 
at the end of inflation. In that sense we will consider mostly the superhorizon modes, or 
late-time asymptotics of the modes: 



p < a{T)H. 



(78) 
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The evolution of the modes at superhorizon scales taking into account the bulk interac- 
tion through the integral kernel is nonlocal. Consequently the evolution (and the possible 
dissipation) of superhorizon modes depends on their past evolution at subhorizon scales 
through the interaction of the subhorizon modes with the bulk. At superhorizon scales the 
bare frequency of the bound mode is imaginary 

^lir) = -Ih\t). (79) 

The bare bound mode coefficient c[ (r) excluding the interaction with the bulk is a super- 
position of a growing mode and a decreasing mode 

^ A+e+t^" + A_e-i^", (80) 

at superhorizon scales and in the WKB approximation. Since we have rescaled the fields 
in ( !53l l by the factor a'^/^(r), the true bound mode is actually a superposition of a constant 
mode and a decreasing mode. The bare solutions ( 1801 ) correspond also to the superhorizon 
asymptotics of the four dimensional solution given by the Hankel function in the limit of a 
pure de Sitter four dimensional geometry, 

— > ~ e "(^)^ (81) 

where rj is the conformal proper time and ai^r) = e^'^ . The decreasing mode is unobservable 
and has little cosmological interest for the formation of the large scale structures in the 
universe. In the following we consider only the dominant solution composed of the growing 
mode. 

We look now for the growing mode of the dressed bound state coefficient with the 
interaction with the bulk taken into account by using the superhorizon Ansatz 



Cb{T) = exp[+7r] (82) 



where 



3 
7 = -H + A-f, A7<if. (83) 

Since the bulk interaction kernel is nonlocal, we should take care when inserting this su- 
perhorizon Ansatz into the integral kernel in eqn. (l75ll because the rigorous complete 
integration should take into account the contribution of the subhorizon regime, where the 
mode has a form different from the superhorizon Anstaz ( l82l) . However, we conjecture 
that, when observing the dressed superhorizon solution today, the integrated effects of the 
subhorizon modes in the brane-bulk interaction are subdominant compared to the contri- 
bution of the superhorizon modes. This is because the amplitude of subhorizon modes is 
suppressed compared to the amplitude of superhorizon modes. The choice of the Ansatz 
( l82l) at any time for solving the integro-differential equation (l75l ) for the bound state is 
equivalent to smooth the weak amplitude subhorizon oscillations of the complete mode 
solution. In that sense we consider the superhorizon Ansatz ([82l) sufficient to solve the 
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integro-differential equation and compute the order of magnitude of the dissipation rate of 
the bound state at superhorizon scales. 

We now proceed to compute the attenuation factor of the growing bound mode in 
the regime Hd ^1. In this regime the curvature radius of the extra dimension is much 
larger than the Hubble radius so that the dissipation rate of the bound state into the extra 

the 



dimension is expected to reach its maximal value. As we discussed in section 16.11 
plateau potential is a suitable approximation in this regime for the study of the brane-bulk 
interaction. The height of the plateau potential is then ¥^{7)/^ ~ AH'^ij) and the dressed 
bound mode coefficient evolves at superhorizon scales according to the equation 

9 f^°° 

h{T) - -H\T)Cb{T) + / dsK{s)Cb{T -S)=0, 

^ Jo 

where K{s) is the bulk interaction kernel taken in the superhorizon limit, eqn. 
Inserting the Ansatz ( l82i ) into ( l84l) yields the algebraic equation 



(84) 



Zl. 



7 



9 

4^ 



H' 



128 



H^H 



n 



2^,2 



dk- 



k'-f 



(P + AHY {k^ 



im 



y) 



which to linear order in A7 becomes 
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such that 
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The minus sign in A7 induces attenuation of the growing mode. The plateau potential 
approximation has the advantage of eliminating the processes of diffraction or reflection in 
the bulk by smoothing the bulk inhomogeneity, keeping only the dissipative processes due 
to the extra dimension. Even if the bound state did not dissipate at subhorizon scales and 
had only its frequency shifted, the bound mode would be observed to dissipate classically 
in the extra dimension at superhorizon scales as 



Cb(T 
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where c^ (t) is the bare growing mode, 
state at superhorizon scales is given by 



r(r) 



The dissipation rate r(r) 
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(89) 



The emitted bulk gravitons oscillates at all times because their frequency remains above the 
plateau of the potential. We may also express the attenuation factor in terms of the four 
dimensional slow-roll inflation parameter en = —H/H'^, characterizing the adiabaticity of 
the expansion of the universe, and the number A^ of e-folds as 



-jr{r)dr 
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(90) 



The graviton bound state decays nonlocally and quadratically in the slow-roll factor. 
The dissipation of the zero mode of the graviton during inflation is thus subdominant 
compared to the dissipation of the inflaton perturbation (section [4]) . 
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7 Scalar perturbations in the approximate background: 
adiabatic perfect fluid on the brane 

In this section we study the evolution of scalar metric perturbations in the approximate 
background geometry flSOll accurate near the brane. These perturbations couple to the 
adiabatic perturbations of a perfect fluid on the brane. They are all described by the 
single master field Q as previously discussed in section [3l We take the Fourier transform 
in the three spatial transverse directions because of the homogeneity and the isotropy, and 
separately evolve each Fourier mode. In the approximate background geometry (iSOl) near 
the brane, the equation of motion (l26l) of the master field in the interval < .^ < +00 is: 
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where p is the transverse momentum. Similarly to the tensor case in section EJ we reduce 
the equation of motion f lOTl) . rescaling 
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to the equation 
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which contains terms dependent on ^ representing the inhomogeneity of the AdS bulk. The 
equation of motion is not separable yet. However, as in section [6] for the tensor case when 
we introduced the plateau potential, we again simpliplify the dependence on ^ retaining 
only the leading behavior of the coefficients in ^ about ^ = 0. We set 
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reducing f l93ll to 
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where 
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The bulk curvature i is contained in A and the intrinsic brane curvature H is contained in 
uj. We highlight the "wrong" sign in the potential A(r) for the scalar perturbations since 
H{t) < 0. The brane behaves as a repulsive potential for the scalar perturbations so that 
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there is no scalar gravitational bound state near the brane as long as there is no coupling 
to matter on the brane. 

The case of a perfect fluid on the brane (with equation of state Pm = wpAi) gives rise 
to some difficulties in computing the evolution of perturbations because the expansion of 
the brane is no longer adiabatic. H/H"^ is not small except when the fluid is subdominant 
compared to the brane tension a (cosmological constant). Here we explore the damping of 
the fluid perturbation on subhorizon scales because in this regime the adiabatic approxi- 
mation of the expansion is valid. At subhorizon scales the geometry of the brane appears 
quasi-Minkowskian. 

The boundary conditions (l27l) linking the master field to the fiuid perturbations on the 
brane become in the approximate background geometry ([50l) near the brane 
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It is possible to simplify the boundary conditions ([971) for the scalar metric perturbations 
by expressing them in terms of the gauge-independent matter perturbation {p — a) A = 
6p — 3H6q. One obtains the following boundary condition 



n' 



<y2{T) 



Q 



— K 



:(P 



a] a" 



5=0 p^ 

which, in terms of the rescaled master field (l92l) . becomes 

.2 (P - ^) «'/' 



A(r). 



(9^ + A(r))fi|^=o = -/^' 



pz 



A(r), 



(98) 



(99) 



where A is given in (I96j) . 

The equation of motion for the matter perturbations on the brane can also be derived 
from the boundary conditions ( l97l) by imposing an equation of state on the matter pertur- 
bations. We consider adiabatic matter perturbations with 
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where the sound speed of the fiuid is c^ = P/p. One obtains the equation of motion for 
the matter perturbation A, similar to the equation obtained in [2T] . 



A+ (8 + 3cl - 6e) HA 



p^cl 



(l0 + 6c2-14e + 3e2) 



2i 



(5 + Scl - 9e) 



12 



A 



(101) 



where according to the braneworld dynamics [4],[6| 
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p = -3H{P + p) = -3H{l + w)pM, 

e = 1 + w, pm = p-a, Pm = P + a. (102) 

The equation of state of the fluid on the brane is Pm = wpM- We may also rescale the brane 
degree of freedom as e^^i'^i jv''^'^''^~^v "^A ^ A to obtain an oscillator equation. Therefore 
the equations coupling the brane scalar perturbations to the bulk scalar perturbations may 
be summarized in the generic form 
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where A and on are given in (J96l) and 
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The immediate observation is the absence of time derivative couplings for the brane-bulk 
system in (11031) . Consequently the dissipation of the brane fluid into the bulk can only 
be nonlocal despite the perfect fluid being local {i.e. truly localized on the brane). Say 
otherwise, no local friction term such as r(r)9^ can appear in the effective equation on 
the brane, this means that the dissipation of the brane fluid contains a memory or a time 
delay. The nonlocal brane-bulk interaction comes from the curvature effects. 

The underlying physics can be illustrated using the following mechanical model. Con- 
sider a string coupled to a harmonic oscillator at the boundary with two springs as indi- 
cated in Fig. El The equations of this model are similar to those obtained for the scalar 




q(t) C) 




Figure 6: Mechanical model for braneworld scalar perturbations. The string, representing 
the bulk gravitons, is coupled by an intermediate spring to a harmonic oscillator representing the 
brane fluid. 
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braneworld perturbations 

il - u" + ujl^i^{t)u = 0, 
7{t)u%=o - uJc(t)u\x=Q = -ujl{t)q{t), 

q+{ujl{t)+ujl{t))q = ul{t)ul=o, (105) 

where 7 is the string tension expressed as a frequency, and cUq and ul are the spring 
constants of the two springs expressed as frequencies squared. A toy model very similar to 
this mechanical system was considered in [27|. The intermediate spring between the mass 
and the string plays the role of a buffer. In the absence of the buffer the coupling with the 
string would bring a first derivative local term (friction term) in the effective equation for 
the harmonic oscillator, indicating a local dissipation of the oscillator down the string. 
With no uniform restoring force on the string, the equations would be 

ii(x, t) — u"{x,t) = 0, 
-ju^x = 0,t) — ujc'^u{x = 0,t) = —uj^qit), 
q{t) + {uJo^ + uj,'')q{t)-uj,\{x = 0,t) = f{t) (106) 

where f(t) is a forcing term. For coefficients not depending on time, the ansatze u(x,t) = 
u ■ exp[iuj{x — t)], q(t) = q ■ exp[— zcjt], and f(t) = f ■ exp[—iujt] give 

- u^q{uj) + (c^o' + uJc^)q{uj) - , \ q{uo) = f{u;), (107) 

which may be rewritten as 
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and A = oj^ j^. Consequently, eqn. p08l) becomes 
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qit) + (cuo' + ^c)q{t) - Lo/X / ds exp[-As] q{t - s) = f{t). (110) 

In the limiting case where A ^ ujq.Uc, we may approximate 

/"OO 

X ds exp[-As] q{t - s) ^ q{t) - A"^g(t), (111) 

Jo 

so that eqn. (II 101) may be approximated as 

m+im+^o'qit) = fit) (112) 
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We next consider the wave equation for the string including a restoring term, so that 
there is a minimum frequency oOmin for the propagating modes. One has 



u-u + u^i^u = 0. 

In this case the kernel is modified to become 
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When ujmin is not too large, oJmin < A = Ci;^/7, there is an isolated pole at cj = — zA where 
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(115) 



and a branch cut on the real axis extending from oj = —Umin to a; = +ujmin- To obtain a 
causal kernel, we must take a contour passing above the branch cut on the real axis. For 
t > 0, where the kernel has support, we evaluate the contour by deforming it so that there 
are two contributions K(t) = Ki{t) + K2{t) (see Fig. [7](a)). The first term originates from 
the contour encircling the pole in the clockwise sense 



K^(t) = -(cu,'A)e(t)exp[-At]^^' ^'^"^'" 
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The second originates from the contour encircling the branch cut in the clockwise sense as 
indicated in Figl7|(a). 
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Figure 7: Analyticity properties of the kernel in the complex plane. Panel (a): if 
i^min < ^ = ^c/l then there is an isolated pole on the lower imaginary axis and a branch cut 
[—comin, -'r^min] on the real axis. Panel (b) and (c): if LOmin > A = co'^/j then the pole moves to 
the real axis under the branch cut as well as its mirror image coming from the second Riemann 
sheet. We may deform the branch cut by pushing it to negative infinity (panel (c)). 

As Umin approaches A = uj'^/'j, the pole on the lower half-plane approaches u = 0. 
There is a second pole at the mirror image a; = +iA. However, this pole lies on the other 
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Riemann sheet hidden beyond the branch cut. When Umin exceeds A = cij^/7, these two 
poles collide and scatter at right angles so that they lie on the real line on the branch cut 
(or just below after the deformation in Fig. [7](b)). We can make both these poles disappear 
by deforming the branch cut downward as indicated in fig. [7](c). The branch cut has been 
pushed all the way to negative infinity. Now the kernel integral may be expressed as 
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where the contours C3 and C4 are indicated in[7](c). We may rewrite this as 
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In fact this expression also holds for ujmin < A = uj'^/'y. At large times the kernel ( 1119^ 
behaves asymptotically as 

r^m *^~ ^n/T\ , cos (cUmmt) - sin (Wmint) ,^^^. 

K{t) ^ C(l)7Vw^m -^^ ^ -■ (120) 

We also compute the kernel in Fourier space. With the forcing term /(t) = /■exp[— zcjt], 
the oscillator equation is 

ijj ^A 
- w2g(a;) + (^0^ + ^/)g(^) - - — . .% ^ =g(^) = /(w), 

A - i^^^ - ui^^ 

G~\u)q{uj) = fico). (121) 

The time average of the power transmitted to the oscillator (which is a measure of the 
dissipation) is 

= u mpiq] - ^miq]) . 

= \f\^uQ \g{u)\ . (122) 

If the kernel has no imaginary part, the system exchanges energy but does not dissipate 
since the average power vanishes. For uj < Umin the imaginary part of the kernel vanishes. 
From G{uj) we deduce the average dissipated power 
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A2 (coi - co^f + (a;2 - ulJ H + u;^ _ ^2f 

The dissipated power is zero when oj"^ < ojf^^^. One can understand qualitatively how the 
behavior changes as uimin is varied. The string may be regarded as a high pass filter with 
threshold oOmin- Only radiation of higher frequency can escape down the string to infinity, 
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thus leading to a net flux of power away from the oscillator. At lower frequencies the 
string is excited around the oscillator but in a localized way with the amplitude decaying 
exponentially with distance from the oscillator since the extra momentum k = \/up- — co'^^^ 
from the bulk wave equation becomes imaginary. In this case there is no dissipation in the 
long term because the string vibrates in phase with the oscillator. Of course the coupling 
shifts the frequency of the oscillator because the string adds both inertia and restoring force 
to the oscillator. 

For the brane-bulk system, the threshold is given by 

^Ln = ^'+^' (124) 

where A, u have been computed in (l96l) . For a weak coupling to the bulk {H <^ i/^), the 
frequency of the brane degree of freedom is only slightly displaced from its bare frequency 
uq given in eqn. (1104^ . At subhorizon scales {p ^ aH), 

0<uj',^'^ + o(h) < ^Ln^^^,+0(H) (125) 

because c^ < 1. Therefore the dissipated power is exactly zero. In the adiabatic approxima- 
tion, valid on subhorizon scales, the fluid cannot dissipate into the bulk because the natural 
frequency of the bulk mode at the same wavenumber is too high to have a resonance. 

8 Discussion and conclusions 

In this paper we presented analytic approaches to estimate the order of magnitude of 
the dissipative effects in braneworld cosmology affecting the evolution of both scalar and 
tensor perturbations on the brane. Because of the non-uniform expansion of the brane, 
the degrees of freedom on the brane interact with the bulk gravitons, possibly leading to 
dissipation from the four dimensional perspective. These effects are encoded in the AdS 
bulk retarded propagator. From the dressed brane propagator, obtained by resumming 
the bulk backreaction effects at all order in the brane-bulk coupling, we may extract the 
effective dissipation rate of certain brane degrees of freedom. We assumed an adiabatic 
expansion on the brane (if ^ if ^) in order to apply our methods. 

Our analytic results provide more intuition on the previous numerical simulations re- 
garding braneworld slow-roll inflation [23]. For the scalar perturbations with a slow-roll 
inflaton field on the brane we obtained, without any other approximation, the local dis- 
sipation rate of the inflaton perturbation due to its interaction with the bulk gravitons, 

r ~ if (if/if^ j (if£/Vl + W^P^, which is linear in the slow-roll factor. This is our main 

result and agrees with the numerical results obtained by Koyama & al. in [23]. In the 
high-energy limit of braneworld inflation (if^ ^ 1) the correction to standard inflation due 
to the coupling to bulk metric perturbations is found to be linear in the slow-roll factor as 
well. Moreover the dependence on Hd of the dissipation rate correctly fits the behaviour 
of the slow-roll correction plotted in [23] . 

For the tensor perturbations we encounter purely nonlocal brane-bulk interaction, that 
is why we simplified the backgound geometry: we explicitely neglected the nonlocal pro- 
cesses on the brane due to backscattering of bulk gravitons in curved AdS assuming that 
the dissipation is the dominant effect on the brane by considering the near-brane limit 
of the bulk geometry (plateau potential). This approximation is well legitimized in the 
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high-energy regime of braneworld inflation {Hi ^ 1) because the bulk curvature becomes 
negligeable compared to the acceleration of the brane. We were able in this way to compute 
the slow-roll correction to tensor perturbations due to the coupling of the zero mode with 
the continuum of bulk gravitons. We found a small attenuation of the growing mode of the 
graviton bound state at superhorizon scales and at high energy due to its interaction with 
the bulk modes and we obtained that this correction is quadratic in the slow-roll factor, 
r ~ H{H/H'^Y. This is a new result. The inflaton field perturbation thus decays at a 
larger rate T ~ H{H/H^), linear in the adiabatic slow-roll parameter, than the decay rate 
the graviton bound state T ~ H{H/H'^)'^, quadratic in the slow-roll parameter at high- 
energy. This difference is the consequence of the locality of the interaction between the 
inflaton and the bulk modes as opposed to the nonlocal interaction between the graviton 
bound mode and the bulk modes. We also found that an adiabatic perfect fluid on the 
brane does not dissipate into the extra dimension at subhorizon scales because the mini- 
mum frequency of the bulk modes at the same wavenumber as the fluid is too high to have 
a resonance. We encounter some difficulties in applying analytic methods for the perfect 
fiuid at superhorizon scales because at such scales the expansion of the brane is no longer 
adiabatic, at least when the fiuid dominates the tension of the brane. 

It would be interesting to go beyond the plateau potential approximation by taking into 
account the inhomogeneity of the AdS bulk and the diffraction of the gravitons by the bulk 
which lead to nonlocal effects on the brane. The nonlocality is encoded in the AdS bulk 
propagator and depends strongly on the bulk curvature. In that sense the exact expression 
of the AdS propagator is necessary to compute the nonlocal effects. The curved geometry 
of the bulk also requires the existence of metastable gravitons or quasi-bound states on the 
brane even if the brane is static as it was found by Seahra [2H|- We could compare the fiux 
of the tunneling bulk gravitons with the fiux of the quasi-bound states. This will be the 
object of a future publication. 
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